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Abstract 

In this letter, a Hamiltonian of the impurity model is constructed within 
the framework of the open boundary Heisenberg XYZ spin chain. This 
impurity model is an exactly solved one and it degenerates to the inte- 
grable XX Z impurity model under the triangular limit. This approach 
is the first time to add the integrable impurities to the completely 
anisotropic Heisenberg spin model with the open boundary conditions. 

As quantum-mechanical model of magnetism, Heisenberg's model [jlj is very fruit- 
ful in the theory of magnetism and is actively studied after Bethe's work on isotropic 
case of the XXX model. Yang and YangjpJ] generalized Bethe's method to the 
XX Z model. Baxter in his remarkable papers gave a solution for the com- 
pletely anisotropic XYZ model with the use of the Bethe ansatz method. Faddeev 
and Takhtajan proposed the quantum inverse scattering method for the XYZ spin 
model and simplified Baxter's formulae [Q. Then, many exactly solved models have 
been proposed and solved by the coordinate Bethe ansatz ||, the functional Bethe 
ansatz || [7] and algebraic Bethe ansatz method || |J, etc. Recently the greatest 
progress has been made for the quantum impurity problems such as Kondo problem 
and tunneling in quantum wires for the one dimensional electron systems. We know 
that the impurities play an important role in the strongly correlated electron sys- 
tems and even a small amount of defects may change the properties of the electron 
systems. Then it is very important to construct the integrable systems including 
the impurities. The pioneering work on the impurity model with the integrability 
was carried out by Andrei and Johannesson [ IT]] for the isotropic Heisenberg chain. 



It were extended to the Babujian- Takhtajan spin chain in Refs. [II]. Bediifig 



E/3Zer and Frahm[T2] solved the integrable model with the impurity coupled with 



periodic t — J chain fll"3|, |I~4"|, |1"5|. Schlottmann and Zvyagin have introduced the im- 



purity in supersymmetric t — J model via its scattering matrix with the itinerant 
electrons [T7j . The Hamiltonian of the system and other conserved currents 



can be constructed in principle by the transfer matrix. They have discussed also 
the magnetic impurities embedded in the Hubbard model [T^] and a finite concen- 
tration of magnetic impurities embedded in one-dimensional lattice via scattering 
matrices [[BJ. 

Exactly solved systems with the open boundary conditions have been studied 
earlier in Refs. [20-24] and a general approach to construct open quantum spin 
model is given by Cherednik[25] and Sklyanin|26[]. Then many exactly solved model 
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with the boundary conditions have been proposed and solved after these pioneering 
work [27-42]. We notice that the impurities may cut the one-dimensional system 
into the small part when they are introduced and (then) the open boundary systems 
are formed with the impurities at the ends of the systems. Then the integrable 
impurity model [43-49] can be constructed from the open boundary system. In this 
paper, we devote to construct an integrable Hamiltonian of the impurity model for 
the completely anisotropic XY Z Heisenberg spin chain where the impurities are 
coupled to the ends of the system. As is well known, the XY Z spin model is 
described by the Hamiltonian 



N-l 
n=l 



(1) 



where J±, J 2 and J 3 are constants; the spin operators a ] n have the form 

a{ = I <g> • • • <g> a j ® • • • <g> / (j = 1, 2, 3; n = 1, • • • , N) 

where o" J at site n are the Pauli operators. This Hamiltonian denotes a one- 
dimensional quantum mechanical model of the ferromagnetism. There are N spins 
labelled by n = 1,2, •••,7V on a line. Every spin is associated with the three- 
dimensional vector 7? = (a 1 , a 2 , a 3 ) of the Pauli matrices. The interactions between 
the neighboring spins are expressed by the constants Ji,2,3- The R matrix of the 
system has the form 



it: (A) 



/csnr/snAsn(A + rj) \ 

snA sm/ 

sn?7 snA 

sn(A + 77) 



(2) 



sn(A + if) 



V fcsnr/snAsn(A 
which satisfies the Yang-Baxter equation 

R12 (Ai — A 2 ) R13 (Ai — A 3 ) R23 (A 2 — A 3 ) 

= R23 (A 2 — A 3 ) -R13 (Ai — A 3 ) R12 (Ai — A 2 ) 
and has the properties 

#12 (A) #12 (-A) =p(A)J, 

R 12 (0) = sn^P 12 , R 1 ^ (A) = R 12 (A) , 

R\\ (A) R\\ (-A — rf)—p(X + r))I 

with p (A) =sn 2 r]— sn 2 A. Another kind of the basic properties in the above Yang- 
Baxter equation is the difference property that the R matrices rely on only the 
differences of the corresponding spectrum parameters, which is useful to construct 



(3) 



(4) 



2 



the models with integrable impurities. The boundary K matrices satisfy the reflec- 



tion equation [25, 26 



R12 (Ai - A 2 ) K- (Ai) R w (Ai + A 2 ) K- (A 2 ) 

K- (A 2 ) i? 12 (Ax + A 2 ) ir_ (Ai) R 12 (Ai - A 2 ) (5) 



#12 (-Ai + A 2 ) if+ (Ax) i? 12 (-Ax - A 2 - 2?7) K*+ (A 2 ) 
= (A 2 ) i?i 2 (-Ai - A 2 - 2 V ) kl (Ai) fl 12 (-A x + A 2 ) (6) 

1 2 

where K±= K± ® idv^ an d K±= idy x ® K±. Here I follow the same notations as in 
Refs. [26-29]. The transfer matrix t (A) is defined as 

t (A) = tr {K + (A) T (A) K_ (A) T _1 (-A)} (7) 

with the use of the boundary K± matrices and the monodromy matrix T (A). It has 
the commuting property 

[t(A),t(A')] = 0, 

which can be obtained from the Yang-Baxter equation and the boundary reflection 
equations ( ||) and (||). This property ensures the integrability of the system with 
the open boundary conditions. For convenience, we use the notation 

4 

R(\)=Y J Wj (A) a 3 <g> a 3 (8) 

3=1 

where 



Wi (A) + W 3 (A) = sn(A + 77), 

Wa (A) - W 3 (A) = snA, 

W l {\) + W 2 {\) = sn V , (9) 

W^i (A) - W 2 (A) = fcsn7/snAsn(A + 77). 

In order to add the magnetic impurities to the XYZ spin chain, we define that 

T a (A) = R aR (A + c R ) R aN (A) R aN -i (A) ■ ■ ■ R al (A) . (10) 

Then the inverse of the monodromy matrix is 

T a 1 (~ A) = tt r-R a i (A) R a2 (A) • • • -R a AT (A) i? a fl (A — Cr) . 

P N (A) p (A - en) 
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We choice that the boundary K matrices as 



tf_(A) = R aL (A + c L ) R aL (A - c L ) , 
X + (A) = 1, 



(11) 



which satisfy the reflection equations. Here the impurity parameters cr and cl are 
introduced. They are related to the exchange constants between the particles of 
the system and the impurities situated at the ends of the chain. These exchange 
couplings between the boundary spins can be turned away from their bulk values by 
variation of the parameters cr and cl respectively. It occurs since the construction 
is based completely on the difference properties of the solutions of the Yang-Baxter 
equation (3). The transfer matrix is 



t (A) = Tr {T a (A) K_ (A) T' 1 (-A)j . 



Then we have that 



Tr a 



dR aR (A + c R ) 



RaN (0) RaN-1 (0) • • • R al (0) K_ (0) T" 1 (0) 



A=0 



2p(c L ) * dWj (A + c R ) 

p(cr) h~ 



W 3 (-c R ) , 



(12) 



A=0 



Tr a R aR (c R ) 



dR a N (A) 



d\ 



RaN-i(0)---Rai (0)K_(0)T~ L (0) 



A=0 



P(cl 



snr/p 



^— - {2JiA R a l N a l R + 2J 2 B R a 2 N a R + 2 J 3 C R cr^cr| 
\ C R) 



+ 2J 3 Y,W 3 (cr)W,(-cr)\, 
j'=i 



where 



(13) 



(C fl ) W 3 ("Cfl) + ^3 (C R ) W 2 (-Cr) , 



(14) 



Br = W 1 (cr)W 3 (-Cr) + W 3 ( C r)W 1 (-Cr) 

+W 2 (c R ) W a (-cr) + W 4 (cr) W 2 (-cr) , (15) 

Cr = W 1 (cr)W 2 (-Cr)+W 2 ( C r)W 1 (-Cr) 

+ W 3 (Cr) W 4 (-Cr) + W 4 (Cr) W 3 (-Cr) . (16) 
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By the use of 



P, 



(A) 



nn+l 



d\ 



1 + ksn 2 r] 1 1 1 — A;sn 2 77 



2 2 



A=0 



we have that 



cnrjdnr] 3 3 cn^dm? 
H o ^n+i H o ' 



JV-1 r 

Tr a R aR (c R ) R aN (0) • • • R aj+1 (0) 

7 = 1 I 



dR a j (A) 



dA 



A=0 



•^i(0)---i? a i (0)K_ (O)T- 1 (0)} 

2p(c L ) ^ 
snr] j . 



iV — 1 

n+1 



where J\ — (1 + ksn 2 rj) /2, J 2 = (1 — ksn 2 rj) /2 and J 3 =cnr/dn?7/2 ; and 

1 



(A) p (A - Cfl ) 

p / n x (A + C R ) 



Tr a {T a (0)K_(0)R al (0) 



d\ 



X=0J 

dW 3 (A - c fl ) 



dA 



A=0 



where 



Tr a T a (0) 



(A) 



rfA 



t- 1 (0) 



A=0 



= 2A L <j\o\ + 25^^ + 1C L o\o\ + D L 



2_2 



Ar. = 



dW 1 (A + c L ) 



d\ 

dW 2 (A + c L 



W 4 (-c L ) + 



dW 4 (A + c L ) 



A=0 



+ 



d\ 


A 


dW x (A - 


cl) 


d\ 


A 


dW 2 (A - 




dA 



^3 ("Cl) 



dA 

dW 3 (A + c L 



W 1 (-c L ) 



A=0 



W 4 (C L ) + 



^3 (C L ) 



rfA 

dW 4 (A - c L ) 



dA 

rfW 3 (A - c L ) 



W 2 (-c L 

A=0 

(cl) 



A=0 



W 2 (cl) , 



A=0 



Br = 



dW 2 (A + c L ) 



dW A (A + c L ) 



d\ 

dW l (A + c L 



d\ 



W 4 (-c L ) + 

A=0 A=0 

dW 3 (A + c L ) 



+ 



d\ 

dW 2 (A - c L ) 



W 3 (-c L ) - 



A=0 



d\ 

dW 1 (A - c L ) 



W 4 (c L ) + 



d\ 

dW 4 (A - c L ) 



W 2 (-c L ) 
W 1 {-c L ) 



A=0 



d\ 



dW 3 (A + c L ) 



W 3 (c L ) - 



d\ 

dW 3 (A - c L ) 



A=0 

W 2 (c L ) 



A=0 



A=0 



d\ 



A=0 



d\ 

dW l (A + c L 



W 4 (-c L ) + 



rfW 4 (A + c L ) 



A=0 



+ 



rfA 

dW 3 (A - c L ) 



^2 (~c L ) - 



dX A=0 

dW 2 (A + c L ) 



A=0 



rfA 

dW l (A - c L ) 



W4 (C L ) + 



rfA 

<2W 4 (A - c L ) 



m (c L ) , 

W 3 (-c L ) 
W 1 {-c L ) 



A=0 



rfA 

Dl = j2\- w ^ x + Cl) 



W 2 (c L ) - 



d\ 

dW 2 (A - c L ) 



A=0 

W 3 (c L ) 



A=0 



A=0 



d\ 



d\ 

W 3 {-c L ) + W 3 (c L ) 



W 1 (c L ) , 



A=0 

dWj (A - c L ) 



A=0 



d\ 



A=0- 



Therefore we get that 
snr/ (it (A) 



4p (c L ) rfA 



7V-1 



A=0 



n=l 



2p(c L ) L 2p(c R ) L p{c R )^ 
+ (N - 1) J 3 . 



(22) 



(23) 
(24) 



£ ^ ( Cfl ) ^ (- Ci? ) 



Set 



with 



if = 



snr] 



cit (A) 



A=0 



4sn 2 r] — 4sn 2 Ci dA 
/io = A^J 3 + ^sn^sn (2i]) (l — /c 2 sn 2 r7sn 2 c_L^) 



(25) 
(26) 
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1 — k sn (77 + cl) sn (7/ — cl) 

' • .,, ' . ■■ ) • (27) 



ysn 2 ?7 — sii 2 cl sn 2 // — sn 2 CR 
We get the following Hamiltonian of the impurity model, 

N-l 



n=l 
+ - 



- SllVr > 



sn 2 ?7 — sn'CL 



sn 2 ?7 — sii^cr 



- — 2 — {^i "^ "/? + Jr^nVr + Jr^n^r} > ( 28 ) 



where 



1 + ksn 2 rj 1 — ksn 2 r] cni]dni] 

Ji= 2 ' 1/2 = 2 ' 1/3 = ^ _ ' 

J Lj i = isn 2 r]cnc L dnc L [1 + ksn (77 + c L ) sn (77 - c L )\ , 
■h,2 = ^sn 2 r]cnc L dnc L [1 - ksn (77 + c L ) sn (77 - c L )\ , 



Jl,3 = -sn?7sn (277) (1 — A; 2 sn 2 ?7sn 2 Ci 
• 1 — k 2 sn 2 (77 + Cl) sn 2 (77 — cl 



And the exchange constants Jr^ have the same expressions as the constants Jla ( 
i — 1, 2, 3 ) except for the substitution of the parameters by cr. It is an integrable 
impurity model. The two impurities are coupled to both of the ends of the completely 
anisotropic Heisenberg spin chain. When we take the triangular limit k — > ( the 
supplementary modulus k' tends to 1 ), the elliptic functions sn77,cm/,dmi become 
sin u, cos u, and 1, respectively. The above Hamiltonian reduces to 

1 f^ -1 

H = 2 1 51 ( a n^n+l + ^n+l + COS T]<y^ l <J^ l _^^ 



, n=l 

sin 2 7/ cos c L fix. 2 2 , cos 77 3 3 



sin 2 77 — sin 2 cl V " cos 

sin" 77 — sin" cr \ " " " " coscr 



Sin 7] COS C R f j 1 ,22, COST? 3 3 \ . 9Q N 

m 77 — sin cr V cos Co / 



It is just the Hamiltonian of the impurity model related to XX Z spin model j|7|, ffl^ . 

As the conclusion, we have added the magnetic impurities to the edges of the 
completely anisotropic Heisenberg XYZ spin model. The model is exactly solvable 
and the Hamiltonian of the impurity model is obtained explicitly. The interactions 
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between the impurities and the electrons are described by the arbitrary parameters 
cl and cr, which comes to the difference properties of the spectrum parameters of 
the R matrices of the spin chain. Under the triangular limit, our Hamiltonian of 
the impurity model reduces to the one related to the XX Z spin mo del ^]. We 
known that the XYZ spin model is equivalent to the eight vertex model || and 
some progress has been made for the eight vertex models with the open boundary 



conditions pi], p2[ . Batchelor et aZ|4(| have obtained the surface free energy with 
the boundary K matrix []53|. It is an interesting subject to discuss the impurity 
effects in the two dimensional lattice model in statistical mechanics. Furthermore, 
the Bethe ansatz equations of this system with the impurities can be obtained by 
the use of the general Bethe ansatz method. We may study the thermodynamics of 
the XYZ model including the magnetic impurities and the impurity effects in the 
ground state and excited state. All of these are in investigations. 
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